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Abstract 

We report determination of branching fractions for the decays ip(2S) — > h + J/ip, where h = any, 
7r + 7r~, 7r°7T , 77, 7T°, and 77 through Xc0,i,2- These measurements use 27M tp(2S) decays produced 
in e + e~ collision data collected with the CLEO detector. The resulting branching fractions and 
ratios thereof improve upon previously achieved precision in all cases, and in combination with 
other measurements permit determination of B(x c j 7«W0 an d B(ip(2S) — > light hadrons). 
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The study of charmonium has entered an era in which many of the broad features are 
finally known, increasingly focusing attention upon the details of both production and de- 
cay [1]. The decays of if)(2S) in particular have become very well studied. About 4/5 of 
all ip(2S) decays are through de-excitation, mostly by hadronic transition to the J/ip, but 
also through radiative decay to the Xcj states. Measurements of rates for such processes 
enable meaningful comparison with theory and extrapolation to similar mechanisms in the 
T system. The listed transitions can be used to isolate and study lower-lying cc states. 
For example, the decay ip{2S) — > 7r + 7r~ J/ip is useful because this is the most common and 
the most accessible experimentally; inclusive and exclusive J/ip decays are often identified 
by tagging the recoil dipion. The corresponding ip(2S) decay rate is therefore a target for 
continuing precision improvements. Similarly, Xcj mesons can often be tagged for study by 
the transition photons. The inclusive rate for transitions to the J / ip is also a crucial input 
to predictions for cc annihilation rates because it limits the remainder. Hence refining the 
rates for ip(2S)-to- J/ip transitions remains an experimental priority. The first comprehen- 
sive look at all transitions simultaneously by a single experiment was reported by CLEO 
in 2005; absolute measurements were limited by a 3% uncertainty in the number of ip(2S) 
produced, and many ratios of rates were limited by statistics. Further investigations with a 
larger dataset and improved systematic uncertainties are certainly warranted. 

This article describes a measurement of ratios of branching fractions ip{2S) — > h + J/ip, 
where h = 7r + 7r~, 7r°7r°, r](— > 77, 7r + 7r~7r°), tt°, 77 through 7Xcj (for which process we 
will also use the expression 7(7-//'0)xcj); an d the inclusive branching fraction B{ip{2S) — > 
any + J/ip), in which only the J /if) is identified. We separately determine £>+_ = B(if>(2S) — > 
n + n^J/ip) with the J/ip decaying inclusively, which we denote by J/ip — > X. 

We use e + e~ collision data at and below the ip(2S) resonance, E cm = 3.686 GeV (/ Cdt = 
53.8 pb -1 , corresponding to 27M ip(2S) decays) and E cm = 3.670 GeV ("continuum" data, 
J Cdt = 20.6 pb^ 1 ). Events were acquired at the Cornell Electron Storage Ring (CESR) [3] 
with the CLEO detector mostly in the CLEO-c configuration (95%), with the balance 
from CLEO III. The detector features a solid angle coverage of 93% for charged and neutral 
particles. The charged particle tracking system operates in a 1.0 T magnetic field along the 
beam axis and achieves a momentum resolution of ~ 0.6% at momenta of lGeV/c. The 
Csl crystal calorimeter attains photon energy resolutions of 2.2% for E 1 = 1 GeV and 5% 
at 100 MeV. 

For the measurement of branching fractions relative to one another, the J/ip is identified 
through its decay to or e + e~. We require |cos# t rk| < 0.83 for both lepton tracks, 

where the polar angle 9 is measured with respect to the positron direction of incidence. The 
ratios of calorimeter shower energy to track momentum, E/p, for the lepton candidates, 
taken to be the two tracks of highest momentum in the event, must be larger than 0.85 for 
one electron and above 0.5 for the other, or in the case of muons smaller than 0.25 for one 
and below 0.5 for the other. In order to salvage lepton pairs that have radiated photons and 
would hence lose too much energy to remain identifiable as a J/ip, we add bremsstrahlung 
photon candidates found within a cone of 100 mrad to the three- vector of each lepton track 
at the interaction point (IP). The J/ip candidate is retained only if constrained fits to the 
two tracks and bremsstrahlung candidates to a common vertex and to the mass of the J/ip 
fulfil Xv/d.o.f. < 20 and xli/ d.o.f. < 20, respectively. The xivi/d.o.f. restriction corresponds 
roughly to demanding that the dilepton mass lie between 3.03 and 3.16 GeV. 

For ip(2S) — > any + J/ip, cosmic ray background is rejected based on the distance of the 
track impact parameters to the event interaction point (< 2 mm), and on the J/ip momentum 
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(pj/ip > 50MeV/c). To suppress background from continuum reactions we require pj/^p < 
570 MeV. Radiative lepton pair production and radiative returns to the J /if) (e + e~ — > 7 J /if)) 
are suppressed by requiring \cos9j/^\ < 0.98 and for the dielectron mode by demanding 
cos# e + < 0.5. Decays of if)(2S) to final states not involving a J /if) can contaminate the 
any + J /if) mode if two oppositely charged particles satisfy the lepton identification and 
kinematic criteria. Monte Carlo (MC) studies indicate such backgrounds are very small, 
leading to no background subtraction and an assignment of a systematic uncertainty in rate 
of 0.2% Gu + /i") and 0.1% (e + e~). 

For 7T + 7T~J/'if) [7i°ir J/ifj] we demand the presence of two oppositely charged tracks [two 
7T° candidates]. In both cases, the mass recoiling against the dipion is required to be 
3.05-3.15 GeV. To suppress background from radiative Bhabha and dimuon pairs, where 
the photon converts and the resulting tracks are associated with the tt + 7t~, we require 
m(n + n-) > 0.35 MeV. 

Any photon candidate reconstructed in any of the exclusive modes must lie in the central 
region of the calorimeter (|cos# 7 | < 0.75) and be isolated (not aligned with the initial 
momentum of a track within 100 mrad and not closer than 30 cm to a shower that is 
matched to a track). We impose mode-dependent requirements on the photon energies, E^. 

All 7T° candidates are formed using two photon showers, each of at least 30 MeV in energy, 
that together have an invariant mass of 100-160 MeV. For Ti Q n°J/ip and (7i + 7i^7i°) ri J/if), but 
not for ir°J/if>, we constrain their four-momentum to the 7r° mass. For it J /if), the softer of 
the two photons cannot lie within = 100-200 MeV in order to suppress background from 
7X c j decays. 

To select (Tr+Tr^vr )^ events, we require two oppositely-charged tracks and a tc° candidate, 
which together have an invariant mass of 535-560 MeV. For (77), 7 J/V ; > we require that 
both photons have E 1 > 200 MeV and 771(77) = 500-600 MeV. Leakage of the high-rate 
process n iT J/ip into the other modes with only two photons and no tracks is suppressed 
by requiring that for n°J/if> and JXcj the third-largest shower energy in the event does not 
exceed 50 MeV. Both rj decay modes are combined into a common rj measurement using 
their branching fractions 

For 7T°J/if) [r)J/if)], we only keep events where J /if) candidates have a momentum within 
490-570 MeV/c [170-230] MeV/c. 

For 7(7^/V')x c j' we tally yields, cross-feeds, and backgrounds separately for the following 
windows applicable to -E 7 -i ow , the lower photon energy of the two: 90-145 MeV (Xc2), 145- 
200 (xci), 200-245 (between Xci an d Xco), and 245-290 (xco)- For the last two categories, we 
subject the event to a final kinematic fit constraining the J /if) and the two photons to the 
if)(2S) four- momentum. The J /if) momentum must lie in the range 250-500 MeV to eliminate 
ir° J /if and (77) 7? J/'0 cross-feed, and for J = 1,2, the invariant mass recoiling against the 
two photons must lie near the J /if) mass, 3.05-3.13 GeV. 

After these requirements, the event samples are clean. This is demonstrated for any+J/if>, 
7T + 7T~J/if) and 7r°7r°J/if), as well as i]J/if> and Ti°J/ip in Figs. [U [21 and [31 The remaining 
background is readily estimated and subtracted. The most important backgrounds come 
from cross-feed among our signal modes, for which we account using the measured yield from 
our sample as background normalization. We also calculate background from rjJ/if> with r\ 
decaying other than to 77 and 7r + 7r~7r , normalizing with branching fractions determined 
elsewhere jEf. We subtract background from continuum by counting the yield observed in 
our continuum data, scaled by luminosity and the 1/s dependence of the cross-sections. The 
yields are listed in Table [D 
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We determine the detection efficiency for all modes we study using signal Monte Carlo 
(MC) samples generated with Evtgen 0], including photon production in the decay of 
the J/ip ("decay radiation") [7|, and a GEANT-based [8j detector simulation. The irTtJ/ip 
samples were produced using the EvtGen model VVPIPI, with a slight reweighting of 
the dipion mass distribution to better represent the measured spectrum [2j. Allowing for a 
relative D-wave component in the dipion system of the strength found in Ref. 0] changes 
the efficiency by at most 0.1% (relative). The angular distributions for the Xcj decays were 



generated according to the formalism presented in Ref. [lOJ. In addition, for the (77) Xc7 
modes, we simulated the Xcj line shape as Breit-Wigner distributions out to 20/228/19 
(J = 0/1/2) times the full width, usin g m asses and widths from Ref. [HI ]. and scaled by 
i?7-iow as appropriate for El transitions [12j . We compute the efficiency for any + J ftp from 
the weighted sum of the individual signal MC samples, where the weight is the relative 
occurrence as measured in our data. We will show later that the sum of our exclusive 
channels describes the population observed in the inclusive selection adequately. 

The following contributions to the uncertainty were evaluated: MC statistics (at the level 
of 1%), tracking (0.3% per track, added linearly), photon detection (0.4% per photon, added 
linearly), MC modeling (about 1%, additionally 2% for the two-photon recoil restriction for 
the Xcj modes) trigger (0.1-0.4%, depending on decay mode), uncertainty of the background 
subtraction (stemming from all statistical uncertainties involved). 

The Xco background merits some discussion. The signal rate is heavily affected by the 
understanding of the area between the Xci and Xco, since whatever populates this region also 
feeds into the Xco window. We noted in our 2005 analysis Q an excess of unexplained events 
in Ky-iow — 200-245 MeV and confirm a similar production rate here with more data (Fig. H]). 
However, our simulation of this region has improved in several ways that significantly raise 
the MC expectation for -E 7 -i ow = 200-245 MeV, in particular through the aforementioned 
E^_ low weighting and through the continuation of the Breit-Wigner distribution from the 
Xci across this region and extending into that of the Xco- The Xa,2 measurements are not 
affected by these considerations due to their much larger rate (the Xc2,i,o raw yields compare 
roughly as 20:40:1). The events with -E 7 -i ow = 200-245 MeV exhibit the expected behavior 
for 77 J/^ events and are qualitatively comparable to those in neighboring regions in i? 7 _i ow , 
namely: the quality of the J/ip fit is comparable, there is no indication of missing energy 
or momentum, and the lateral shower profile for photon candidates is consistent with this 
assumption. 

We have studied several ways to explain the data in the intermediate region as well as the 
adjacent areas: variation of the Xci width from 0.92 MeV to 1.00 MeV, removal of the 
weighting, allowance for a phase-space-like component where the ip(2S) de-excitation takes 
place through two non-resonant photons (denoted by (77)111-), and combinations of these. 
Most of the scenarios can give a satisfactory description with different component weightings, 
but result in substantial variations of the Xco rate due to the different apportionment of 
the observed yields to signal and background. Our data do not allow us to distinguish 
among the possibilities explored, which results in correspondingly large background-related 
contributions to the systematic uncertainty (x c0 : 10%, Xci'- 1-0%, Xc2- 0.8%). In any of 
these scenarios, more events in the Xco region are attributed to background than in Ref. 0] ■ 
We quote the result obtained with our nominal xa width, with weighting, and with 
a phase-space-like component (77) nr , as displayed in Fig. HI Since any sensitivity to the 
possible additional (77) nr component is restricted to the region between Xco and xa (and 
not unambiguous there), we cannot quantify the magnitude of such a contribution any more 
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precisely than to say that a branching fraction up to 1 x 10~ 3 is compatible with our data. 

Since we measure ratios of yields, many systematic uncertainties cancel, most notably 
those related to lepton species and J /if) fitting. To verify this, we compare ratios measured 
using J /if} — > e + e~ with those determined using J/if> — > ji + \T . All ratios are close to unity; 
for the seven ratios involving ir + 7T^J/ifj we compute \ 2 = 3.0 for six degrees of freedom. 

Table [I] shows results for the branching fraction ratios, after combining the two measure- 
ments from J /if) — > e + e~ and Adding all ratios of exclusive decays to the inclusive 
one, + J/V)/(any + J/ip)\, leads to a sum that agrees with unity within 0.9% or la. 
This favorable comparison is an indication that the contribution from modes not covered in 
this analysis (such as 77 J/^ through r) c (2S), n + n~ it J / if) , etc.) is small. 

We now describe a measurement of the if>(2S ) — » 7r + 7r~J/ip, J /if) — ► X production rate. It 
proceeds identically to the one detailed in Ref . |13| , where it was used as the denominator in 
the determination of B(J/if) — > £ + £~). The presence of the J/if> is inferred through the 7r + 7r~ 
recoil mass spectrum. The pion candidates must consist of two oppositely charged tracks 
that obey loose quality criteria, p t > 150 MeV, |cos6>| < 0.83, and m(7r + 7r~) > 300 MeV. 
The recoil mass spectrum, shown in Fig.0 is fit with a signal shape using the m(7r + 7r~)-recoil 
spectrum from tc + 7t~ (£ + £~) data events, selected with the same criteria for the pions as the 
inclusive J/if> — *> X sample (overlaid in Fig. [5]), and a second-order polynomial background 
to extract the yield of 7r + 7r~(X)j/^, events: 3.851 x 10 6 . The fit has a confidence level of 
98.8%. The efficiency determination is potentially hampered by the lack of knowledge of 
many J/if> branching fractions. In practice, however, it is found to be nearly independent of 
J I if) decay mode, although it can be expected to depend weakly upon the track multiplicity, 



as detailed in Ref. [13j . We therefore use a sample approximating the spectrum in data with 
a weighted sum of signal MC of various multiplicities, including well-measured decays like 
J /if) — ► £ + £~ . The agreement between data and the sum of MC predictions thus obtained 
is good [l3j. The detection efficiency determined from the weighted sum of individual 
efficiencies is found to be 40.16%. Systematic uncertainties stem from uncertainties in the 
efficiency weights and modeling, translating into 0.7%, sensitivity to the fit range in the 
signal yield extraction, 0.2%. The result for the number of ip{2S) — > -K+n^J/if) events 
produced is (9.589 ± 0.020 ± 0.070) x 10 6 . 



The number of if>(2S) is measured in a manner similar to that described in Ref. [14 
Hadronic event candidates are identified by requiring three tracks (N tT ^ > 3) and restrictions 
on the following quantities: summed energies from tracks and showers, -E v is ; amounting to at 
least 0.3E cm , and furthermore for iVtrk = 3, 4 a summed energy deposition in the calorimeter 
of more than 0.15E cra , and either this sum less than 0.75E cm or the highest individual shower 
less than 0.75i?bcam- The ^-component of the event vertex, 2 vtx , must be within 5 cm of the 
beam spot center. 

The continuum contamination consists of e + e~ — > £ + £~ (£ = e, /1, or r), light hadron 
production through e + e~ annihilation or two-photon collisions, radiative returns to the J/if>, 
and J/ip decays from the extended tail of the J /if) Breit-Wigner distribution. It is subtracted 
statistically by scaling the yield of events passing the selection criteria from the off-resonance 
data sample. The scale factor between off- and on-resonance yields is the ratio of luminosities 



(measured using the process e + e~ — > 77 [15j) multiplied by a 1/s dependence for the cross- 
section behavior. We subtract non-collision events statistically using an extrapolation from 
the tail of the z vtx into the signal region. Other effects are negligible. 

In a MC simulation we determine the detection efficiency for if>(2S) decay events to be 
76%. The underlying MC generator settings incorporate current branching fraction deter- 



6 



minations [ll[ for ifj(2S), Xcj, and J/ip decays, and for the remainder employs JETSET 16 
The agreement between data and MC simulation can be judged from Fig. [61 

To evaluate the accuracy of the estimated efficiency and background subtractions, we 
explore three different scenarios where we vary the requirements on track multiplicity (> 1 
to > 4) with appropriate background suppression criteria modifying the aforementioned 
energy balances, with detection efficiencies ranging between 91% and 65%, respectively 
We find a variation of 2% relative to the nominal setting, which we take as a systematic 
uncertainty. The largest contributors are the dependence on the trigger requirements (1.8%), 
followed by the tracking and energy settings (0.9%). 

The procedure, when applied to a portion of the earlier-taken CLEO-c ip(2S) data, results 
in a slight reduction in the result compared to the method described in Ref. [13]. This is 
understood to be due to updated settings of the MC generator (thereby modifying the 
detection efficiency), a change in the continuum background subtraction (using CLEO-c 
continuum data), and other improvements in the detector description. 

The number of '0(25') decays thus determined is (27.36 ± 0.57) x 10 6 , with a relative 
systematic uncertainty of 2%. The statistical uncertainty is negligible. This leads to a 
branching fraction B(ip(2S) — > tt + tt~J/iI)) = (35.04 ± 0.07 ± 0.77)%. This result is more 
precise than the one presented in Ref. [2j and also somewhat higher. This is chiefly due to 
an improved background and signal efficiency treatment in the inclusive ip{2S) count. 

The last column of Table [J shows the absolute branching fractions, obtained by multi- 
plying the entries in the penultimate column by £>+_ and appropriate cancellation of cor- 
related systematic uncertainties. One independent set of numbers from the data displayed 
in the table is £>+_ together with the ratios to £> + _. Some derived quantities may be com- 
puted, with the correlations properly taken into account. We find B(n° 'J /ip) /B(r]J /ip) = 
(3.88 ± 0.23 ± 0.05)%. Using B(ip(2S) -»■ iXcj) = (9.3 ± 0.4)%, (8.8 ± 0.4)%, and 
(8.1 ± 0.4)% [lH for J = 0, 1, 2, we obtain B{x c j -> jJ/iff) = (1-35 ± 0.07 ± 0.14 ± 0.06)%, 
(40.5±0.3±1.4±1.8)%, (24.1±0.2±0.9±1.2)%, respectively, where the first uncertainty is 
statistical, the second systematic from this analysis, and the third from the input branching 
fraction. Our result for B(x c o jJ/^A is considerably smaller than found in Ref. 0] and 
consistent with other determinations [ll|. The reduction in this rate is a direct consequence 
of improved treatment of Xco backgrounds. 

With B(if>(2S) -> !Xcj,lVc) and B{ip{2S) -> E + £~) {£ = e,/i,r), our results imply [Tl| 
B(ip(2S) — > light hadrons) = (15.4 ± 1.5)%, or 2.9 standard deviations higher than an 
extrapolation arrived at by scaling B(J/ip — > light hadrons) by the ratio of leptonic branching 
ratios, (12.45 ± 0.35)%. Here, all ip{2S) results are taken from the branching fraction fit 
values in Ref. fll| . 

In summary, we have studied the exclusive decays ip(2S) — > J/ip + h {h = TT + n^, n°n , 
i], 7T°) and ip{2S) — > ^Xcj —> 11 J W transitions, J/ip — ■> e + e _ and fi + fi~, with a similar 
strategy applied to all channels. The analysis is complemented by a study of the inclusive 
mode ip(2S) — > any + J /ip. We have determined branching ratios between exclusive modes 
on the one hand and between exclusive modes and (any + J / ip) on the other. We have also 
measured the branching fraction ip{2S) — > tt + tt~ J /ip using the dipion recoil mass spectrum, 
which facilitates transformation of the ratios relative to tt + tt~ J /if) into absolute branching 
fractions. Further quantities are derived. The precision of all quantities given here improves 
upon previous measurements. The results presented here supersede those from Ref. [2|. 
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TABLE I: For each channel: the number of events observed in J/ip — > /-t + n~ after background 
subtraction and the detection efficiency ratio = e(tp(2S) —> h + J/tp^f 1 )/e(ip(2S) — > any + 
J/ipv+v )| the same for J/ip — > e + e~; the ratio of branching fractions B(\p(2S) —> h + J/ip and 
B(ip(2S) — > any + J/Y>); the same with respect to £> + _; absolute branching fractions. 



Channels 


N>> 




N e 


rt 


B/B any (%) 


(%) 


B (%) 


7T + 7T~ J/l/j 


302030 


0.80 


263372 


1.01 


56.04 ±0.09 ±0.62 


= 100 


35.04 ±0.07 ±0.77 




32249 


0.17 


28746 


0.22 


28.29 ±0.12 ±0.56 


50.47 ± 0.22 ± 1.02 


17.69 ±0.08 ±0.53 




9819 


0.27 


8590 


0.33 


5.49 ±0.06 ±0.09 


9.79 ±0.10 ±0.15 


3.43 ±0.04 ±0.09 


7T°J/V> 


289 


0.19 


238 


0.25 


0.213 ±0.012 ±0.003 


0.380 ± 0.022 ±0.005 


0.133 ±0.008 ±0.003 


7(7^/V')xco 


308 


0.22 


253 


0.28 


0.201 ±0.011 ±0.021 


0.358 ±0.020 ±0.037 


0.125 ±0.007 ±0.013 


7(7^/^x01 


13244 


0.34 


11619 


0.44 


5.70 ±0.04 ±0.15 


10.17 ±0.07 ±0.27 


3.56 ±0.03 ±0.12 


7(7^/V')xc2 


6616 


0.31 


5768 


0.40 


3.12 ±0.03 ±0.09 


5.56 ±0.05 ±0.16 


1.95 ±0.02 ±0.07 


any + J/ip 


676889 


= 1 


466153 


= 1 


= 100 


178.4 ±0.3 ±2.0 


62.54 ± 0.16 ± 1.55 



9 




JA|/ Momentum (GeV/c) 



FIG. 1: Distributions relevant to any + J/ip. Top: polar angle of the positive lepton. For 
e + e~ only, we demand cos# e + < 0.5 to suppress Bhabha events with initial/final state radia- 
tion or bremsstrahlung in detector material. Bottom: J/tp momentum. Solid circles show the 
on-tp(2S) data, dashed histogram the continuum data (scaled by luminosity and 1/s) taken at 
E C ra =3.67 GeV, the solid histogram represents the sum of all MC exclusive channels (scaled to 
match the data in signal modes). Arrows appear at nominal selection values. 
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FIG. 2: Plots relevant to the decay ip(2S) -> vr+vr-J/^ (top) and 4>(2S) -> TT°TT°J/ip (bottom). 
The left plots show the dipion recoil mass spectrum and the right plots the dipion mass spectrum. 
The J/tp candidates in the continuum sample arise from the tail of the ip{2S). Symbols are as in 
Fig. 1. The dotted line represents the simulated signal. 
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FIG. 3: ip(2S) -> r]J/xp, <q -> 77 (top) and 7/ -> vr+vr-vr (middle), and ip{2S) -> tt J/> (bottom): 
The J/ , momentum (left) and invariant mass of the decay products (right). Symbols are as in 
Fig. 1. The dotted line represents the simulated signal. 
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FIG. 4: For the decays tp(2S) — > 77J/1/), S 7 -i ow (top) before applying the full-event kinematic fit 
to the area above -E 7 -i ow = 200 MeV, the di-photon recoil mass for Xc2,i (middle [lower] left for 
J = 2 [J = 1]), and -E 7 -i ow after the full-event kinematic fit (lower right). Symbols are as in Fig. 1. 
The dotted line represents the simulated signal. 
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FIG. 5: The dipion recoil mass spectrum for ir + ir~J/vfj, J/ip — > X. Top: data points (black) 
overlaid with the fit result (solid black curve) obtained using a (scaled) signal shape from ir + ir~J/ip, 
J ftp — > £ + £~, and a second-order polynomial background shape (red dashed curve). Bottom: the 
fractional difference between the fit and the data. 
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FIG. 6: Scaled energy in the event E v [ s (left) and number of tracks (right), for data (dashed black 
line / points) and simulation (solid red line). The right distribution has been obtained with the 
selection based on jV t rk > 1- 
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